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Abstract. Atomic ions are mostly neutralized by small grains (or PAH molecules) in current theories of heating 
and cooling in cool diffuse clouds; in the main they do not recombine with free electrons. This alters the ionization 
balance by depressing n(H’'') and n(He''") while carbon generally remains nearly fully once-ionized: charge exchange 
with atomic oxygen and formation of H 2 and OH also depress n(H^) in partly molecular gas. Seemingly restrictive 
empirical limits on (^h are relaxed and higher values for ((h are favored in a wide range of circumstances, when 
grain neutralization is recognized. Maintenance of the proton density at levels needed to reproduce observations of 
HD requires Ch ^ 2 x but such models naturally explain the presence of both HD and H^ in relatively 

tenuous H I clouds. In dense gas, a higher ionization rate can account for high observed fractions of atomic 
hydrogen, and recognition of the effects of grain neutralization can resolve a major paradox in the formation of 
sulfur-bearing compounds. 
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1. Introduction. 


Several authors (Wolfire et al 


Weingartner fc Praine T ^OOlb ; Welty fc Hobbj , 



have commented recently on a seeming side-effect of cur¬ 
rent theories in which the heating of the diffuse ISM is 
dominated by the photoelectric effect on small grains or 


PAH molecules ( 

Omont, 

1986 


D’Hendecourt & Leger 

1987 

; Lepp & Dalgarno, 

1988b; 

Bakes & Tielens, 

1994 

Weingartner & Draine, 

2001c 

). The usual processes of di- 


rect radiative and dielectronic recombination of atomic 
ions become secondary to neutralization during gas-grain 
interactions. When gas is somewhat denser or shielded 
from the interstellar diffuse soft X-ray flux, the proton 
density is much lower than would otherwise be the case 
and the total ionization fraction is much more nearly that 
given by the amount of free carbon in the gas alone; car¬ 
bon remains mostly once ionized. 


Thi s neu tr alization process w a s discu ssed by [Draine fc 
Sutiiil (1987]), Lepp & Dalgarno (1988a) and jLepp et ^ 


(1988) and is not limited to diffuse gas or hydrogen ions. 
Indeed, ion stage ratios like Mg II/Mg I, C II/C I etc. 
should be entirely controlled in this way if large molecules 
cum small grains heat the neutral gas. All this notwith¬ 
standing, recognition of the importance of grain neutral¬ 
ization has occurred only with great hesitation, obser- 
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vational demonstrations of its relevance are still largely 
lacking and it is typically noted by observers (if at all) 
as something of a curiosity: Welty fc Hobbs] ( ^OOlj ) is an 
obvious exception. The original authors were concerned 
with establishing the existence of the effect and determin¬ 
ing the abundance of the responsible agents, postulating 
what they described as ad hoc modifications of the chem¬ 
istry. Now, with the benefit of hindsight and with general 
acceptance of the heating effect of the same population of 
small grains, it is appropriate to recognize that the higher 
ion-neutralization rates associated with gas-grain interac¬ 
tion have other, largely unexplored consequences for our 
understanding of the interstellar medium (ISM). 


The approach taken here to studying these effects is 
largely based on observations of various forms of hydro¬ 
gen because they are uniquely sensitive to the cosmic-ray 
ionization rate in gas shielded from the diffuse soft X-ray 
flux. Determination of the proton density, cosmic-ray ion¬ 
ization rate and efficacy of grain neutralization are inextri¬ 
cably intertwined. We show that some important historical 
limits on (the cosmic-ray ionization rate of hydrogen, 
defined as the primary rate per neutral atom) must be 
restated, because they were derived only in terms of ex¬ 
plicit ratios Cn/agas, where (symbolically) Ogasj the rate 
for gas-phase recombination with an electron, grossly un¬ 
derestimates the actual rate at which ions are neutralized. 
Relaxation of these limits allows recognition that essen¬ 
tially all observations of the various forms of hydrogen 
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(protons, HD, H 2 and in more diffuse gas; H I in dark 
molecular gas) suggest a much-enhanced cosmic-ray ion¬ 
ization rate. 

The combination of faster neutralization of atomic ions 
- which carry the charge in diffuse gas - and faster cosmic 
ray ionization has strong implications for the chemistry 
of diffuse, partly molecular gas. Unlike atoms, molecular 
ions are dissociated by recombination more rapidly than 
they are neutralized by grains. If the cosmic-ray flux is 
enhanced, producing molecular ions more rapidly, without 
a concommitant increase in either the electron fraction 
or the molecular ion destruction rate, the abundances of 
molecular ions may be much enhanced: HD and H^ may 
coexist in the observed amounts, even in gas of fairly low 
density. 

In the next Section, we lay out some of the underly¬ 
ing ionization and ion neutralization mechanisms in the 
context of recent two-phase heating/cooling models of the 
ISM. In Sect. 3 we discuss some received limits on the pro¬ 
ton density and cosmic ray ionization rate in atomic gas, 
and demonstrate that limits on ((h from radio recombina¬ 
tion lines are not restrictive if grain neutralization occurs; 
we conclude that ionizing fluxes could be much higher 
than previously assumed. In Sect. 4 we explore some con¬ 
sequences of the combination of grain neutralization and 
higher Ch for the chemistry of diffuse clouds and we show 

as 


that these include a much-enhanced abundance of Hj~, 


has recently been s een in the ISM ( McCall et ah , 2002 ; 
Geballe et ah , 1999| ), in the same gas which harbors HD. 
Section 5 is a brief discussion of processes in dark, fully 
molecular gas. 


2. Physical models of ionization, recombination 
and the ISM 

2.1. Ionization in two-phase models of heating and 
cooling in the ISM 

Our discussion is based on a calculation of heating and 


cooling very similar to that set out by Wolfire et al 


a spirited defense of classical notions of two-phase equilib¬ 
rium, in the presence of other influences, has recently been 
given by Wolfire et al.| ( ^0031) . This model has been us ed by 
us to discuss H 2 and CO formation in diffuse clouds ([Liszt 


fe Lucas, 2000), the spin temperature of H I in warm gas 


(it is not thermalized by collisions, see jLisztj (|2001|)) and 


the formation of H 2 in damped Lyman-a systems ( [Liszt 
200 ^) . 


Such models embody three primary sources of ioniza¬ 
tion and heating; soft X-rays, optical/uv photons and cos¬ 
mic rays. The X-ray flux is harsh but fragile and easily at¬ 
tenuated by relatively small gas columns. For the reference 
model a ” bare” intrinsic spectrum is attenuated by a nom¬ 
inal gas column (following Wolfire et al. ( 1995| ); Nw 
« 10^® cm“^) whose value can be adjusted to reproduce 
the observed local mean density of electrons in intercloud 
gas (which we will also reference as pertaining to ’Tree- 
space”), (ne) « 0.02 cm“^ (Taylor & Cordeg, 1993). It is 


Table 1. Primary cosmic-ray ionization rates^ and unde¬ 
pleted relative abundances 


Species 

X 

C 

(10-i"s-i) 

[X/H] 

H 

1.00 

1.00 

D 

1.00 

1.5 X 10"® 

He 

1.08 

7.9 X 10"^ 

H 2 

2.00 


C 

3.83 

3.6 X lO"'* 

N 

4.50 

1.0 X lO"'* 

0 

5.67 

7.5 X lO"'* 


_ from _ the _ ^UMIST reaction database 

(http://www.rate99.co.uk) 


a property of two-phase models that the electron density 
in free space varies little with the total hydrogen density 
n(H) for n(H) < (ne)/^c where fc is the abundance of free 
carbon in the gas phase. 

The ionization rate of hydrogen required to ionize the 
intercloud gas at the observed levels is 10“^^s“^, about 
100 times larger than the presently accepted ionization 
rate due to cosmic rays. This intercloud ionization rate 
has changed little over many years ( Field et al. , 1969 ) and 
is probably the best-determined upper limit on acceptable 
Ch in unshielded regions. In cases with substantial inter¬ 
vening gas columns, the X-ray flux is propagated accord¬ 
ingly in our models; the spectrum hardens as the overall 
flux attenuates. 

The photoionization rate can be parametrized in terms 
of a scaling parameter GO, relative to a reference interstel¬ 
lar radiation field, but this will not be varied here. When 
gas parcels are modelled, the attenuating column densities 
are averaged appropriately over the assumed gas geome- 
t ry. For H 2 and C O, appropriate shielding factors are used 
( Lee et al. , 1996 ). In cool diffuse gas, the equilibrium tem¬ 
perature and pressure at a given density are quite sensitive 
to GO, and to the assumed depletion of coolants (mainly 
carbon). 

The reference value for the primary cosmic ray ion¬ 
ization rate of hydrogen atoms is here taken as Ch = 
10“^^s“^. The ionization rate of a hydrogen molecule is 
twice this and the total primary destruction rate of an 
H 2 molecule due to cosmic rays, including dissociation, is 
2 X 1.08 X Ch- Table 1 gives the cosmic-ray ionization rates 
of various species scaled to t he reference value of Ch , t aken 
from the UMIST database { http://www.rate99.co.uk ), as 
are all the rate constants cited here, unless otherwise 
noted. 


2.2. Grain neutralization 

The small grain population acquires electrons from the 
ambient gas by charge exchange with (neutralization of) 
positive ions having a superior ionization potential (typ¬ 
ically, 6.8 eV according to Omont (1986)) and loses elec¬ 
trons through the photoemission processes responsible for 
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heating. In this way, the grain heating rate is dependent 
on the ionization level, and both are controlled by the 
incident ionizing fluxes and particle density. 

The neutralization rates follow from the grain charg¬ 
ing formalism of praine fc Sutin (1987), integr ating over 
the grai n size spec trum, assume d here , as in Bakes fc 
Tielen^ (|l994| ) and Wolfire et al. (1995 ) to be an exten¬ 


sion of the classical MRN spectrum ( Mathis et al. ^ 1977 b 


see the l ast pa ragraph in this subsection. [Weingartner fc 
Draine| (|2001a| ) give analytic approximations to the low- 
temperatu re cha rging rates (not used here) and ( Lepp fc 
Dalgar R 11988 ^) give a rate constant for cold gas. Grain 
neutralization rate constants for various ion species dif¬ 
fer only inasmuch as their thermal velocities differ, i.e. all 
are proportional to the inverse square root of the atomic 
weight. 

The abundance of small grains needed to heat the gas 
is about 100 times below the Solar abundance of carbon. 
Thus, in diffuse gas where the ionization fraction is at least 
equal to the relative abundance of free carbon (typically 
40% of Solar) and grains lose electrons rapidly to photo¬ 
electric emission, the charge is carried by atomic ions and 
free electrons. In dense molecular gas where photoemis¬ 
sion is nil and the equilibrium gas-phase ionization frac¬ 
tion falls to values below 10“^, the negative charge on 
grains plays a substantial role in the overall charge bal¬ 
ance ( Lepp fc Dalgarno , 1988a ). 


Weingartner & Draine (2001c) have recently revisited 
the grain spectrum and its effects on heating and cooling. 
They derive smaller sticking coefficients for a somewhat 
more copious small grain population and quote heating 
rates which are typically 20%-50% larger than those used 
here. By contrast, Wolfire et al. (2003) recently reduced 
the grain abundance and heating rate by a factor two in 
their reference model for local gas. In light of this conflict¬ 
ing guidance, we simply note that our assumptions appear 
to be intermediate. 


2 . 3 . Charge exchange between oxygen and hydrogen 


The most rapid gas-phase process affecting protons in 
atomic gas is not recombination with electrons (or grain 
neutralization) but charge exchange with neutral oxygen. 
As originally described by Field fc Steigman ( 1971 ) in the 
context of a purely atomic gas, the effect of this charge ex¬ 
change was_to_chaiigethe_Oj||_^(Oj|atio^ 2 mMOji^ITj|J^(n[^ 


I) to e3:p(-232 K/T) n(H+)/n(H I), where the exponential 
represents the slight difference between the ionization po¬ 
tentials of H and O (the latter’s being larger) and the 
higher ionization fraction of oxygen, absent charge ex¬ 
change, follows from the 1971 version of Table 1. Basically, 
oxygen is expected to be ionized much more rapidly than 
hydrogen for a given cosmic-ray flux but is forced to do¬ 
nate much of the ionization to hydrogen. The effect on H 
was presumed to be relatively unimportant, owing to the 
smallness of the 0/H ratio. 



Fig. 1. The ratio of the effective gas-phase recombination 
rate n(e)a^^^ to the total rate at which H+ ions are neu¬ 
tralized by gas-phase recombination and grain neutraliza¬ 
tion, in two-phase equilibrium. The strength of the ion- 
grain interaction is artificially lowered in steps of 2 to 
l/64th of the reference value. Shown are results for Solar 
and depleted ((ic.N.o = 0.4) gas-phase abundances. 


Perhaps for this reason, O 


not explicitly mentioned by Wolfire et al. 


H charge excha nge was 
(it was 


included). However, the displacement of charge is accom¬ 
panied by a slight diminution of n(p) in our models and 
inclusion of charge exchange becomes important when the 
proton density must be determined, especially if the abun¬ 
dance of molecular hydrogen is considered. In even slightly 
denser gas having a molecular hydrogen component. O’*" 
ions interact rapidly with H 2 to form OH;^ and this, like 
all molecular ions, recombines much more rapidly with 
electrons than does any once-ionized atom. In the absence 
of charge exchange with hydrogen, the equilibrium abun¬ 
dance of 0+ would be several orders of magnitude lower 
in gas having a molecular hydrogen fraction of even a few 
percent. 

As noted by O’Donnell fc Watson ( 1974 ) and Barsuhn 


& Walmsley (1977), maintenance of the O"*' fraction can 


actually siphon charge from hydrogen when typical diffuse 
cloud chemistry operates more strongly in denser but still 
diffuse gas. This is another effect of molecules (the first 
being grain neutralization) which tends to lower the pro¬ 
ton abundance in diffuse gas. Although it is possible that 
a proton charge transferred to oxygen is returned to the 
gas after some chemical chains, this will not generally be 
the case. The upper limit on the neutralization rate of pro¬ 
tons due to oxygen charge exchange is the rate at which 
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0+ interacts with H 2 to form 0H+ (the rate constant is 
1.7 X 10“®cm^s“^) and our models simply assume that the 
process is either off (in so-called atomic gas) or on at this 
rate whenever the abundance of H 2 has been calculated. 
As shown below, converting 25% of a moderate density 
cool gas to H 2 results in a decline of 30%-50% in n(p). 


Stancil et al. (1999) give the rate constant for the 
exothermic O -I- H charge exchange channel as 3 x 
10“^°cm^s“^below lOOK, increasing as for T > 100 
K. Increasing this rate has no strong direct effect on n(p) 
- charge exchange already dominates - but would increase 
n(0''") somewhat in more strongly molecular regions, re¬ 
moving slightly more charge from hydrogen. 


2.4. Particulars and parameter sensitivities 


The reference model conside red here assumes the s mall 
grain properties described by Bakes fc Tielenj ( 1994 ); we 
do not adjust the grain properties except artificially to 
demonstrate their specific influence on the gas. Figure 1 
illustrates the efficacy of H-ion neutralization by grains, 
as used in the published calc u lation s of tw o -phase heat - 
ing and cooling Wolfire et al. ( 1995| ); Liszt ( 2001 , 20021 ). 
In this hgure the per-grain cross-section for neutraliza¬ 
tion is artificially varied away from the reference value 
(labelled 1); other aspects of the influence of small grains 
are not disturbed {e.g., the heating rate). The calculation 
is for free space, with and without gas-phase depletion of 
coolants (0,C) and charge sources (C). Grain effects be¬ 
come stronger as the gas-phase ionization level declines 
(as could be the case in more heavily shielded regions). 
Clearly, the usual gas-phase recombination processes play 
a distinctly secondary role in cooler gas and the total rate 
at which protons are neutralized is much higher than has 
typically been assumed previously. 

Shown in Fig. 2 are the electron and proton densities as 
functions of the total density, when the grain neutraliza¬ 
tion cross-section is altered between I and 1 /32 times the 
reference value, as in Fig. 1. In the top panels, no elemen¬ 
tal gas-phase depletion is inclu ded (corresponding most 
nearly to the phase diagrams of Wolhre et al. ( 1995 )); be¬ 
low, it is taken as a typical value. In the right-hand panels, 
the rate of primary cosmic-ray ionization of atomic hydro¬ 
gen is increased by a factor of 16 from its reference value of 
1X 10“^^s“^ (all cosmic-ray rates scale in the same propor¬ 
tion; the rate of cosmic-ray ionization of O I, for example, 
is some 6 times that of H I, see Table 1). The straight line 
in each panel represents n(C“'')= 3.6 x 10“^(5c 11 (H) and 
departures of n(e) from it are due mostly to n(p) (He“'" may 
contribute somewhat disproportionately at lower density, 
owing to its high cross section for ionization by X-rays). 
In the lowest two panels, further complications are intro¬ 
duced; the gas is shielded from soft X-rays by half the 


hydrogen column density of a typical H I cloud (Spitzer, 
196S, 1978) and, for n(H) > 100 cm“^, a moderate frac¬ 


tion of the gas (one-fourth of the H-nuclei) is assumed to 
be molecular. 



10 100 , 10 100 
n(H) [cm-3] 


Fig. 2. Electron and proton density in two-phase equilib¬ 
rium for nominal (Ch= 10“^^s“^; left) and 16 times en¬ 
hanced (right) cosmic-ray ionization rates under a variety 
of conditions. At top, for free space with Solar abundances; 
middle, free space with typical ISM depletions; bottom, 
with depletion, a substantial intervening neutral column 
density (Nw) shielding the soft X-ray flux and, at high 
density (n(H) > 100 cm“^), one-fourth of the H-nuclei in 
H 2 . In each panel the strength of the grain neutralization 
interaction is artihcially scaled from 1 times to 1 /32 times 
the reference value (steps of 1/2). In each panel the pro¬ 
ton density n(p) (shown red in appropriate media) is the 
family of curves which lies lower at right. 


In cooler gas, neutralization of protons by small grains 
causes n(e) to decrease by factors of 2-3 under reference 
conditions,and n(p) decreases somewhat more. The hydro¬ 
gen emission measure product n(e)n(p) typically declines 
by a factor 10-20 in cool material, as suggested by Fig. 
1. The proton density follows Ch closely only when X- 
ray ionization is weak; otherwise, n(p) increases only very 
sluggishly with increasing Ch. 
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3. Ionization balance in atomic gas 

3.1. Limits on Ch in cool gas from radio recombination 
lines 


Radio recombination lines (RRL) of hydrogen offer the 
possibility of fairly direct observation of the proton abun¬ 
dance. To sense the cosmic ray ionization rate, it is neces¬ 
sary to avoid lines of sight having any contribution from 
H II regions and to select gas which is sufficiently shielded 
that the X-ray ionization is weak. The electron density 
has a large contribution from n(p) at all densities in un¬ 
shielded regions if ion-grain neutralization is neglected, 
but only at low densities when it is included. In practice 
it has turned out to be difficult to avoid gas having a large 
molecular fraction, although the complications introduced 
in this way are never taken into account by observers. 

Extant limits on the cosmic-ray ionization rate in cool 
atomic gas arise from (many) unsuccessful attempts to de¬ 
tect low-frequency (decametric) radio recombination lines 


(RRL) of atomic hydrogen (Shaver, 1976; Payne et al 


1984, 


1994). The basic analysis begins with the condition 


of detailed balance for hydrogen 


n{e)n{p)ai^'^ = — n{p)) 

where Th is the total ionization rate and is the 
usual gas phase radiative recombination rate into the n 
> 2 levels, so that no ionizing photons are inadvertantly 
discounted. The strength of a RRL is proportional to the 
line of sight integral f n{e)n{p)dl = < n{e)n{p) >L (L 

is some length), and to a generalized function of the prop¬ 
erties of the host medium and the atom, accounting inter 
alia for non-LTE conditions among the level populations 
in the recombined atoms. The strength of the 21cm line 
is assumed to be proportional to <n(H I)>L = <(n(H)- 
n(p)>L over the same path and to a much simpler func¬ 
tion of temperature (assuming that the spin temperature 
is thermalized, a good assumption in cooler gas). In this 
case the equation of detailed balance is symbolically recast 
as 


RRL strength Th 

21cm strength 


F(ji{H), n(e), T, i/, etc.) 


The function F may have zeroes for certain combina¬ 
tions of density, temperature and quantum level, so the 
RRL to be studied must be chosen with some care. For 
his model of the recombination, Shaver ( 1976| ) derived 
limits Ch ^ 2 — 4 X 10“^®s“^ in a dozen or so H I fea¬ 
tures seen at low latitudes in the inner Galaxy. Many of 
these are fairly distant (as judged by their velocities) and 
most are sources of CO emission in galactic surveys, so 
they must be fairly dense; neglect of molecular processes 
somewhat overstates the sensitivity of the RRL measure¬ 


ments to Ch- Payne et al. (1984) set more stringent limits 
Ch < 0.4 — 1.6 X 10“^'^s“^ in gas toward 3C123 which is 


also a well-studied molecular line source (prutcher , 1980D . 


These limits have been considered a significant demon¬ 
stration that the cosmic ray ionization rate is incapable 
of ionizing the intercloud gas to the observed degree (see 
Sect. 2.1). 

In the context of the current models for ion neutral¬ 
ization, such limits must be revised upward to the same 
extent that the effective recombination rate increases (Fig. 
1), which in general makes them uninteresting. Including 
grain neutralization reduces n(p) by factors of 6 or so at 
n(H) ^ 30 cm“^, and n(e) by factors of 2-3, which is suffi¬ 
cient to ensure that the sought-after hydrogen RRL could 
never be observed. The RRL data really cannot be used to 
derive interesting limits on Ch, if it is accepted that grain 
neutralization is important. 


Eventually it was recognized, largely on the basis of 
the behaviour of low-frequency carbon RRL (which were 
seen), that the function F(n(H),n(e),T, v, etc.) is subject 
to considerable uncertainty in terms of the model which 
is assumed for detailed balance in the recombining atoms. 
Payne et al. (1994) analyzed low-frequency carbon RRL 
and upper limits on the optical depth of the hydrogen RRL 
toward Gas A, and did a quite detailed heating and cool¬ 
ing calculation of two-phase equilibrium (excluding grain 
neutralization processes). Their limit on Ch for one fea¬ 
ture with an inferred density n(H) = 200-500 cm“^, Ch 
< 6 X 10“^^s“^, is perhaps less interesting now than the 
somewhat more robust inference that n(p)/n(e) < 0.3. In 
unshielded atomic gas, such low proton fractions cannot be 
achieved unless grain neutralization dominates. However, 
the gas seen toward Gas A is dense and dark enough that 
i t cannot be considered to be either diffu se, or atomic 
( Wilson et al. , 1993 ; Liszt fc Lucas , 1999| ). The cosmic- 
ray ionization rate Ch is not meaningfully constrained by 
the Gas A data either, if grain neutralization is recognized. 


3.2. Determination of n(e) via optical spectroscopy 


In their work 
tralization, Lepp et al 


introducing the notion 
(119881) 


of grain neu- 
discussed its influence 


on commonly-occuring measurements of ion ratios (Mg 
II/Mg I and the like); the basic idea was that lowering such 
ratios with grain processes mimics the effect of higher elec¬ 
tron and total particle densities. This in turn may help to 
reconcile discrepancies between densities measured from 
ion ratios and those gleaned from study of the excitation 
of G I fine-structure or CO rotational levels. This mecha¬ 
nism has, however, only been accepted very grudgingly by 
observers of ion stage ratios, who are more likely to pay 
lip service than to calculate the effect, and most likely to 
sweep the matter under the rug. 


Ratios of atoms in different stages of ionization depend 
on the total electron density, and would be sensitive to 
large n(e)-n(C+) ~ n(p); this dependence is complicated 
and changed by grain neutralization, of course. If it were 
found that n(e) >> n(C+) in sufficiently dense gas (the 
ionization is supported by hydrogen at low density any¬ 
way), this would suggest that n(p)/n(C+) was large, and 
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n(e) 

0.1 

0.01 


n(e) 

0.1 

0.01 


Fig. 3. Electron density (right scale; red in appropriate 
media) and ion stage ratios n(X II)/n(X I) for X = mag¬ 
nesium (top) and carbon, with substantial gas-phase de¬ 
pletion of carbon and oxygen Jc = (5o = 0.4. The ion stage 
ratios (shown blue in appropriate media) are the families 
of curves which are lower at right; families are created 
by scaling the grain neutralization as in Figs. 1 and 2. 
The ion ratios observed toward C Oph and HD 192639 
(see Sect. 3.2) are indicated by horizontal lines and the 
inferred hydrogen and electron densities are given by the 
vertical line segments connecting the families of curves at 
these densities. 


that Ch was large enough to overcome whatever recombi¬ 
nation/neutralization processes were occurring. In terms 
of Fig. 2, we see that this really cannot happen if grain 
neutralization occurs substantially as indicated, so high 
proton densities in sufficiently dense gas could be a strong 
contraindication of grain neutralization. 

Toward C Oph, Ravage et al. ( 19921) measured N(Mg 


II)/N(Mg I) = 105 and N(Fe II)/N(Fe I) = 537 for the 
denser gas at -15 kms“^, and found n(e) = 0.065 cm“^ 
and = 0.046 cm“^ from the two ratios, respectively, based 
on a rather simple analysis of ionization equilibrium in 
an atomic gas in free space, neglecting grain neutraliza¬ 


tion and cosmic-ray ionization. For n(p) = 0, these values 
of n(e) would be reproduced for n(H) = 300-400 cm“^, 
n(C''')/n(H) = 1.5 x 10“"^. These densities are very nomi¬ 
nally appropriate to the feature in question, perhaps mak¬ 
ing it appear that n(p)/n(C+) is not larg e. How e ver, a 
self-consistent application of the analysis in Spitzer (1978) 
shows that n(p)/n(C“'") would be of order unity even for 
Ch= 10“^^ s“^ and n(H) a factor of two smaller, neglect¬ 
ing grain neutralization. In any case, savage et al. (1992) 
remarked that their values of n(e) were an order of magni¬ 
tude smaller than in the classic analysis by Morton (1975), 
which was based in part on the N(C H)/N(C I) ratio (now 
measured at 60 but previously (ibid) was 30). 

The close agreement in n(e) from Mg and Fe ioniza¬ 
tion ratios was illusory, as relevant f-values for the Mg 
H lines were subse quently revised upward by a factor 2.4 
( Fitzpatrick , 1997| ), lowering the derived Mg H/Mg I ion 
ratio and requiring a substantial upward scaling of the 
electron density derived from it. Ion stage ratios seen to¬ 
ward C Oph are given in Fig. 3 here, which shows n(X 
H)/n(X I) for X = magnesium and carbon in a free-space, 
two-phase heating and cooling calculation assuming a sub¬ 
stantial depletion of gas-phase carbon and oxygen 6 = 0.4 
so that [C]/[H] = 1.44 x 10“"^. 

In each panel of the figure the ion ratios are the fami¬ 
lies of curves which lie lower at right and higher at left; the 
other family in each panel is the electron density. Families 
of curves were created by scaling the grain neutralization 
rates (as in Figs. 1 and 2) from 1 to 1/32 times the refer¬ 
ence value. The shaded trapezoidal outlines in each figure 
illustrate where the models reproduce measured ion ratios 
(the horizontal base) and what electron densities are found 
with and (largely) without grain neutralization. The n(e) 
values cited by Ravage et al. (1992) are reproduced nearly 
exactly for the ion ratios they quoted. 

If grain neutralization is ignored, the hydrogen (elec¬ 
tron) densities corresponding to the Mg and C ion ratios 
are 530 (0.15) and 1140 (0.25) cm“^ respectively. At the 
temperature inferred for this gas from the ratio of popula¬ 
tion in the J=1 and J=0 levels of H 2 , T = 56 K, the ther¬ 
mal pressure found from these ion ratios would be much 
higher (10-20 times) than in ferred from excitatio n of the 
fine-structure levels of C I (Jenkins et ah, 1983) or the 


rotat ional excitation of CO ( Liszt] , i979| ; Liszt fc Lucas 
[1994 ). Including grain neutralization greatly reduces the 
densities, to 105 (0.023) and 145 (0.028) cm“^ , so that 
n(C+) ^ 0.66 n(e) and the n(H)-T products derived from 
different analyses largely agree . 

bonnentrucker et al. ( 2002 ) present new ion stage ra¬ 
tios for several elements toward the low (mean) density 
but translucent {J^b-v = 0.64 mag) line of sight to HD 
192639. Their ratios for carbon and magnesium are re¬ 
peated in Fig. 3; in their analysis n(e) « 0.1 cm“^ for 
either species, slightly higher than here. Along this line of 
sight there is a large discrepancy between the local den¬ 
sity derived from excitation of the C I fine structure levels, 
which indicates n(H) = 10-20 cm“^ for T = 100 K, and the 
substantially higher densities n(H) « 120 — 200 cm“^ indi- 
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cated by the ion stage ratios. Sonnentrucker et al.| ( ^002 ) 
ascribed the seemingly high electron fraction to a mixture 
of H I and H II gas along the line of sight but noted that 
their analysis might be modified by grain neutralization. 
Now we see that grain neuturalization does indeed remove 
the discrepancy. 

Weingartner & Draine (2001b) assessed existing ion 
stage ratios in some clouds of low inferred (fro m C I excita¬ 
tion) dens i ty n(H ) = 10-20 cm“^ observed by [Fitzpatrick 
fc Spitzer ( 1997 ) and found that reductions in n(e) and the 
scatter among the n(e) indicated by differing ions could 
be achieved, but that discrepancies between electron and 
particle densities derived from Ca Ill/Ca II and Mg II/Mg 
I persisted. 

Welty & Hobbs ( ^001 ) discussed a large survey of K 
I spectra and, in the absence of K II spectra, derived 
statistical properties of the cloud ensemble based on a 
mean value for the depletion Jk- The median electron 
density n(e) =0.1 cm“^ is typical of the results discussed 
here. These authors discussed grain neutralization at some 
length l and concluded it wao rooponoible for a factor throe 


E 

Q 

X 



11 Itttt 

increa dc in N(K I), reconciling various prooouro ootimatoo. 


Welty & Hobbs (2001) noted Hobbs’ long insistence that 
n(e) « ii(C“'') for clouds seen in Na I, K I, etc, which we 
know now is only possible in diffuse gas if grain neutral¬ 
ization dominates. 


4. Ionization balance in diffuse gas with 
appreciable H 2 

4.1. Observations of HD 

HD does not partake of the same self-shielding which 
is the sine qua non of high H 2 -fractions. Hence the ra¬ 
tio n(HD)/n(H 2 ) might be expected to be very small in 
diffuse gas, well below the inherent [D/H] ratio (locally, 


Fig. 4. HD and H 2 colum n densities. N(HD) valu es are 
f rom Spitzer et al. ( 1974) exce pt for o Per from Snow 
( 1975 ), C Per from ^now ( 1977 ) and ( Oph from Wright 
fc Morton ( 1979 ). N(H 2 ) values are from Ravage et al. 
( 1977j ). The line labelled n(p) = 0.0035 cm“^ is from Eqn. 


[D/H] = 1.5 X 10 ® as quoted by Moos et al. ( 2002| )). 


That this is not so is due to fractionation and charge ex¬ 
change processes which, however uncertainly, nonetheless 
allow an inference of n(p) without explicit reference to 
the processes by which protons are neutralized (Black & 
Dalgarno, 1973; Watson], 197f; Jura , 1974 ; O’Donnell fc 
Watson, 1974; Spitzei, 1978). 


1, evaluated for T = 75 K and Fhd = 5 x 10 ^^s ^ as in 
free space. 


observed quantities and physical constants, with no ex¬ 
plicit dependence on either the density or recombination 
rates [ibid). We have 


N(HD)/N(H2) Thd,, , ,k2 ,,,n(H2)^ ,41, 

- 2)^] exp (-) 

( 1 ) 

where n(H) = n(H I) -I- 2 n(H 2 ) and the photodissoci¬ 
ation rate of HD Thd = Fh,, = 5 x 10~^ ^s~^ under typical 
conditions in free space ( Lee et al. , 1996] ). The derived n(p) 
is nearly independent of the molecular fraction in the gas 
for k 2 /ki = 2.1 and there is no explicit dependence on 
density (only the H 2 fraction). 

The available (Copernicus) observations of HD in the 
nearb y ISM are summariz ed in Fig. 4; N(HD) valu es are 
from Spitzer et al] (|1974) except for o Per from Snow 


In a purep atomic gas ionized by cosmic ra 3 rs the ion - 


Q ppitzer 

Per: 


izationi and recombination rates of H and D atoms wrould 
be very nearly the same (the grain neutralization of D+ 
is slower by a factor \f2) but a strong and slightly en¬ 
dothermic charge exchange with protons H+-I- D -l-AE ^ 
D+-I- H (rate constant ki = 10“® cm^s“^, AE/k = 41 
K) tends to force n(D+)/n(D 1) « n(H+)/n(H 1) exp(-41 
K/T). In the presence of H 2 a rapid and relatively strongly 
exothermic reaction D+ -I-H 2 ^ HD -f H+ forms HD with 
rate constant ^2 = 2.1 x 10“® cm^ s“^. 

If only charge transfer and H 2 fractionation neutralize 
D'^' a relatively compact expression gives n(p) in terms of 


( 1975|), C Per from Snow (1977) and C, Oph from Wright & 
Morton] (1979]). The iltr error range for ( Per is large, 0.5 


rates in this Section are taken from Table 1 of 


(1998) 


Stancil et al 


dex. With the obvious exception of the three datapoints 
outlined (10 Lac, o Per and ( Per) the regression line fit 
to the data, log N(HD) = —5.73 -I- 0.99 log N(H 2 ), has 
very nearly unit slope corresponding to uniform physical 
conditions and N(HD)/N(H 2 ) « 1.0 — 1.2 x 10“®. From 
Eqn. 1 we find n(p) « 0.0035 cm“^ at T = 75 K in free 
space. 

Thus the nominal proton densities inferred from 
N(HD) are small compared to expected values of n(e) for 
all but the most discrepant datapoints, and may be repro¬ 
duced in slightly darker gas lacking X-ray ionization either 
with nominal cosmic ray ionization rates in the absence of 
grain neutralization, or, with it, if Ch is substantially en- 
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hanced (Fig. 2 at bottom). Earlier analyses attempted to 
infer Ch from n(p) and found an enormous range, for in¬ 
stance from 7 X for ^ Per to 10 

( P ’Donnell fc Watson , 1974). 


4.3. and HD in the context of H 2 


formation 


^®s ^ for 10 Lac 


4.2. Observations of 


The detection of in diffuse gas has had some ironic 
twists. Chemical models employed in the original discus¬ 
sion of grain neutralization (Lepp & Dalgarno, 1988a; 
Lepp et ah, 1988) were poisoned by inclusion of the then- 
fashionable assumption of a very small low-temperature 
gas-phase recombination rate for This led the au¬ 
thors to a gross overprediction of N(H 3 ), which they sug¬ 
gested would be observable; their suggestion seems not 
to have been acted upon in a timely manner, thus de¬ 
priving the world of a seeming corroboration of the in¬ 
correct recombination rate. Instead, the recombination 


rate was corrected 

' Amano, 

00 

1 —1 

Larsson et al.. 

1993 

Sundstrom et al. 

, 

994 ), lowering expectations for the 


Figure 5 shows the run of H 2 , HD and in simple dif¬ 


this was considered surprising. With substantial fluctua- 
tio ns, N(Hj~)/EB vft! 8 x 10^"* cm“^/6 mag (see Fig. 14 


of [McCall et aL| (|2002|) ), or N(H;)')/N(H) « 2.3 x 10"® 
for N(H)/Eb_v= 5.8 x 10^^ cm“^. The clear increase in 
column density with line of sight extinction strongly ar¬ 
gues that H;)" is widely-distributed in a relatively common 
cloud component of the ISM, not segregated in clouds of 
too-high density (which would be relatively rare). 

The abundances of HD and H;)" really cannot be con¬ 
sidered except in the broader context of a model for H 2 
itself (see the next subsection), but the chemistry of H;)" 
is quite simple in broad outline. Very nearly every ioniza¬ 
tion of an H 2 molecule leads to formation of Hj^, which 
in diffuse gas is destroyed mainly by recombination with 
an electron. The rate constant for this process is taken as 
03 = 1.2 X cm®(300/T)°-®s“^, which is our interpo¬ 
lation among experimental values which agree marginally 
at 300 K and have somewhat different temperature depen¬ 
dences. H;)" is destroyed to a much smaller extent (typi¬ 
cally, 100 times) by reaction with C and CO (rate constant 
2 X 10“® cm®s“^) or O (rate constant 0.8 x 10“® cm®s“^). 
The CO fracti o n in d iffuse gas is quite small (c/. Fig. 8 of 
Liszt fc Luca^ ( ^OOPj) ), and grain neutralization at a rate 
l/VS that of H+ is also negligible in diffuse gas. 

The rate equation determining the H;^ abundance is 
typically cast as n(H^) n(e) 0:3 « 2 CHn(H 2 ) but the 
implicit dependences of n(H 2 ) and n(e) upon <Ch pre¬ 
vent n(H|^) from increasing linearly with Ch for large Ch 


fuse cloud models such as we have reported on in Liszt 
^ Lucas ( 2000] ) and Liszt ( 2002| ). A spherical clot of gas 
is immersed in the same conditions for which the calcula¬ 
tions of two-phase equilibrium are performed and molec¬ 
ular hydrogen in the uniformly illuminated gas is allowed 
to self-shield. The H 2 abundance is solved iteratively in 64 
radial zones along with the temperature, electron density, 
HD, and Hg abundances etc.. Represented at top in the 
Fig. 5 is something akin to a standard Spitzer H I cloud; at 
bottom is a model with four times higher number and col¬ 
umn densities. At right, the small grains have been treated 
differently at the same hydrogen number and column den¬ 
sity; at top, their neutralization activity (alone) is sup¬ 
pressed by a factor 32. At bottom they are removed from 
the calculation to this degree. The normalization is differ¬ 
ent for H;)" and HD because N(H 2 ) is generally not known 
along lines of sight where N (Hj^) was detected (only N(H) 
is inferred, from the selective extinction using a standard 
gas-grain conversion). 

Figure 5 shows that the observed abundances of HD 
and Hj^ are both reproduced in the models at left, with 
grain neutralization and a much-elevated (factor of 10-20) 
cosmic ray ionization rate. In the absence of grain neu¬ 
tralization, either the typical HD/H 2 ratio can be readily 
reproduced with the reference value of Ch , or H;)" can be 
reproduced with a much (20-100 times) higher value of Ch, 
at the cost of grossly overproducing HD. 

Because the HD abundance is sensitive to the ioniza¬ 
tion level, it is seen to increase toward the cloud edge (the 
effect is only barely resolved); this would seem a natural 
explanation for the high HD/H 2 ratio toward the low H 2 
column density line of sight toward 10 Lac. It is also pos¬ 
sible that HD is actually on the verge of becoming self¬ 
shielding. Observed values of N(HD) are already larger 
than those of N(H 2 ) at the knee of the N(H 2 ) vs. N(H) 
plot (for instance Fig. 8 of Liszt fc Lucas] ( |2000| )). 


The higher ionization rates seemingly demanded by a 


distributed HC component have been noted by (McCall 


as suggested by Eqn. 3 of McCall et al. ( 2002| ). In our 
model we actually solved a small reaction network for 
the self-consistent determination of the abundances of the 
hydrogen-bearing species, including gas-phase formation 
of H 2 via H~. 


jet al. (|2002| ) and van der Tak fc van Dishoeck ( ^OOPj) . The 
latter authors also derived values of Ch from observations 
of H13CO+ in the dense clouds associated with the same 
protostars used as background sources for the Hj^ spectra, 
quoting Ch= 2.6 ± 1.8 x Noting the discrepancy 

in Ch, they ascribed the increase of H^ with Eb-v to 
intervening clouds of density n(H) ^ 10^ cm“®, in which 
carbon was mainly in the form of either C I or CO. We 
would take strong exception to this view, which does not 
account for HD and H;)" in the same gas, but the cosmic 
ray ionization rate in dense gas was subsequently revised 
by Doty et al. ( p002 ), who quote Ch = 5.6x 10“^’^s“^with a 


factor three uncertainty. This does not exclude the higher 
values of Ch preferred here. The models used by these 
authors do not seem to have included grain neutralization. 
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fractional distance into cloud 


Fig. 5. Calculated HD, H 2 and H;^ fractional abundances for three values of the cosmic ray ionization rate Ch= 
10~^^s~ ^, 10~ ^ ^s~^a nd for two uniformly illuminated spherical models. At top is a ‘standard’ Spitzer 

( ^pitzei , 196S, 1978) H I cloud with n(H) = 32 cm“^,N(H) = 4 x 10^° cm“^. At bottom is a model having 4 
times higher number and column density. At right the effect of small grains is diluted in each model, at top by 
simply scaling down the grain neutralization rate, and at bottom by scaling down the quantity of small grains (and 
so, all associated effects). The normalizations of the relative abundances differ because N(H 2 ) is not known along 
most lines of sight with known N(H;^). N(HD)/N(H 2 ) « 1.2 x 10“® along unexceptional lines of sight (Fig. 4) and 
N(H;^)/(5.8 X 10^^ cm~^FiB-v) ~ 2.3 x 10“® in the survey of McCall et al. (2002). 


5. Implications and considerations for dense gas 

5.1. Metal depletions and the problem of sulfur 
chemistry in dense gas 

The Solar elemental abundances of trace metals (Mg, Na, 
Fe, Si etc.)^ 2 — 4 x 10“^, while small compared to that 
of carbon, are much higher than the implied ionization 
fractions in dense dark molecular gas (see just below). 
Thus, even a small ionization fraction in any one of such 
species could dominate the overall charge balance in dark 
clouds. 


It was noted by Oppenheimer & Dalgarno (1974) that 
molecular ions will efficiently transfer charge to such metal 
atoms in dense gas. Because atomic ions recombine so 
much more slowly, in the gas (so the argument went), the 
overall ionization fraction is much increased in this way (if 
grain neutralization is neglected, as has been the case in 
this subject). In turn, because rather low ionization frac¬ 
tions are required for maintenance of high abundances 


and efficient deuteration in dense gas, chemical modellers 
have typically adopted what is described as a ’low-metals’ 
case, in which many metals are assumed to be strongly 
depleted from the gas phase {6 0 . 001). This includes 

sulfur (for example see Nillson et al. ( ^000| )), although it, 
unlike the other species noted just above, does not un¬ 
dergo charge exchange with the most abundant molecular 
ions (H)^, OH)^ and HCO+). 

Such extreme sulfur depletion leads to several contra¬ 
dictions. Unlike, s ay. Mg, Na or Fe, su l fur is undepleted at 
lower extinction ( Savage fc Sembach , 1996 ) and the pat¬ 
terns of gas-phase depletion overall do not explain why it 
should be missing in denser gas. More importantly, remov¬ 
ing so much sulfur from the gas phase makes it impossible 
to form commonly-observed compounds like CS, SO, etc. 
at the observed levels without harboring abundances of 
OH or O 2 which are 100 or more times what is observed, 

( p000| ). The sulfur- 


see Turner (1995) and Nillson et al. 


bearing compounds are formed by reactions of free neutral 
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sulfur atoms with these species and the abundance of the 
sulfur-bearing molecules scale linearly with the amount of 
free sulfur (ibid). 


Alth ough some of the discussion of Lepp fc Dalgarno 
( |l988a| ) is mooted by recognition that dissociative recom¬ 
bination of H;^ is efficient, and although they did not ad¬ 
dress sulfur explicitly, these authors demonstrated that 
grain neutralization renders the issue of metal depletion 
irrelevant to the free electron abundance in dense, dark 
molecular gas. The charge carried by atomic ions is sup¬ 
pressed by grain neutralization; much of the entire nega¬ 
tive charge contribution resides on small grains (see Sect. 
2.2 here), and the free electron density varies little (less 
than a factor two) whether metals are depleted or present. 
In this case, the abundance of gas phase sulfur may safely 
be assumed to be large enough that models with accept¬ 
ably small OH and O 2 abundances - when they emerge - 
will easily produce enough of the sulfur-bearing molecules. 
Additionally the gas phase metals will be largely neutral, 
owing to grain neutralization. This acts to increase the 
abundances of most sulfur-bearing compounds because S+ 
is less chemically productive. Consequences for the abun¬ 
dances of molecules containing other metals largely remain 
to be explored. 


5.2. The atomic gas fraction in dense gas 


In a dark or molecular cloud, very nearly every cosmic ray 
ionization results in formation of an H;^ ion and the ulti- 


mate appearanc e of q = 2-3 hydrogen atoms ( |Cravens~fc 
Dalgarnol , 1978| ). Hydrogen atoms leave the gas phase by 
sticking to (and forming H 2 on) large grains whose den¬ 
sity is assumed proportional to n(H) = n(H2)/2. The rate 
constant for this process is usually denoted by R and the 
equilibrium density of H-atoms in a dark, fully molecular 
gas is n(H I) « q^n/R, independent of n(H 2 ). Typically, R 
« 3 X cm^s“^ is inferred in diffuse gas ( |Jura , 1974; 

Spitzer, 1978| ; |&y et ah , 2002) bf, so n(H I) « 1 cm 


independent of n(H 2 ), for the reference value of Cu¬ 


lt has been recognized since the work of Allen & 
Robinson ( 1976| ) that the fraction of atomic hydrogen in¬ 
side dark clouds, while usually small, is still often much 
larger than expected for steady-state equilibrium between 
H 2 formation at typical rates in cold gas and a cosmic ray 
ionization rate Ch = 10“^^s“^. Observations of CO emis¬ 
sion and H I self-absorption in low-latitude 21 cm pro¬ 
files by Liszt et al. ( 1981| , 1984) suggested n(H I)/n(H 2 ) 
« 0.005 and similarly high values have recently been cited 
by Willacy et al. ( 2002| ) for some local dark clouds. The 
latter authors invoked a process of turbulent diffu s ion to 
account for the H-fraction, while Allen fc Robinsoij ( 1976 ) 
asserted that dark clouds were not old enough for full H I 
- H 2 conversion (also see Rawlings et al. ( 2002| )). 

The higher values of Ch suggested by this work nomi¬ 
nally imply an enhanced atomic hydrogen fraction in fully 
molecular gas, in agreement with observation. However, 
the rate of H 2 formation is really well-known only in dif¬ 


fuse gas, where photoprocesses dominate and grains are 
incompletely dressed. The rate of H 2 formation in dark 
gas has only been interpreted in the context of the diffuse 
cloud data and/or relatively small Ch- Theoretical mod¬ 
els, which might account for changes in the formation rate 
when photodesorption is nil and depletion proceeds, are 
largely lacking. 


5.3. The cosmic-ray heating rate in dense gas 


The st andard calculatio n of heating in dense dark gas is 
that of Goldsmith (^001) who used a (cosmic-ray) heating 


rate 10“^' n(H 2 ) erg cm“^ s“^. Using the mean value of 
the heat deposited per primary event in a weakly ionized, 
f ully-m olecular gas calculated by [Cravens fc Dalgarno 


( |1978D , AQ= 7 ev, we find that this heating rate cor¬ 
responds to Ch = 4.5 X 10~^'^s 


^ some 4.5 times our ref¬ 


erence value; Goldsmith (2001), without attribution, as¬ 
sumed AQ= 20 ev. Thus, increasing our primary ioniza¬ 
tion rate by a factor of 20 would provide heating which is 
only a factor 4.5 greater than that in Goldsmith ( ^001 ). 
The local cooling function varies as if gas-grain 

coupling is neglected and rather more rapidly if it is in¬ 
cluded. A substantial increase in Ch (compared to our ref¬ 
erence value) could be tolerated without inducing major 
changes in the temperature of dense gas; indeed, to heat 
molecular gas to the observed degree apparently requires 
a value of Ch which is at least 5 times larger than our 
reference value. 


6. Summary 

Gas-phase atomic ions are rapidly neutralized by the same 
population of small grains which heats the diffuse inter¬ 
stellar medium. We considered the effects of such pro¬ 
cesses on various forms of hydrogen; protons in atomic 
gas; protons (deuterons), HD and Hg in diffuse gas hav¬ 
ing a molecular component, and H I in dense molecular 
clouds. 

In atomic gas, limits on the cosmic-ray ionization rate 
Ch (the primary rate per H-nucleus; see Table 1) from non¬ 
detection of low-frequency hydrogen recombination lines 
must be revised upward by such large factors that they are 
no longer restrictive. If grain neutralization is recognized, 
limits on Ch from the RRL are higher than the rate needed 
to ionize the intercloud medium in the absence of the soft 
X-ray flux, as in the original theories of two-phase equilib¬ 
rium. Ion stage ratios like N(G H)/N(G I), N(Fe H)/N(Fe 
I) etc. provide accurate electron densities only when grain 
neutralization is fully taken into account. 

HD and the newly-discovered H^ may be shown to 
exist in the observed amounts, even in classical cool H 
I clouds of moderate density, if grain neutralization pro¬ 
ceeds and the cosmic-ray ionization rate has values Ch 
^ 2 X 10“^®s“^. Models having lower values of Ch will 
reproduce HD if grain neutralization is neglected but fall 
grossly short of explaining H;^. 
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In denser gas, a higher cosmic ray ionization rate 
naturally accounts for the (often) surprisingly high H 
I fractions, without invoking mechanisms such as non¬ 
equilibrium H 2 formation, small cloud ages, or turbulent 
diffusion (all of which might still occur). More impor¬ 
tantly, we noted that inclusion of small-grain charging 
in dense gas removes the need for the ’low-metals’ case 
of severe ad hoc sulfur depletion, thus allowing formation 
of sulfur-bearing molecules without grossly overproducing 
OH and O 2 . 

The present work represents no more than the barest 
first step at including grain neutralization and its neces¬ 
sary consequences (like enhanced values of Ch), into our 
understanding of the interstellar medium. 
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